Quenching Dynamics of a quantum XY spin- 1/2 chain in presence of a transverse field 
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We study the quantum dynamics of a one-dimensional spin- 1/2 anisotropic XY model in a trans- 
verse field when the transverse field or the anisotropic interaction is quenched at a slow but uniform 
rate. The two quenching schemes are called transverse and anisotropic quenching respectively. Our 
emphasis in this paper is on the anisotropic quenching scheme and we compare the results with those 
of the other scheme. In the process of anisotropic quenching, the system crosses all the quantum 
critical lines of the phase diagram where the relaxation time diverges. The evolution is non-adiabatic 
in the time interval when the parameters are close to their critical values, and is adiabatic other- 
wise. The density of defects produced due to non-adiabatic transitions is calculated by mapping 
the many-particle system to an equivalent Landau-Zener problem and is generally found to vary as 
l/y/r, where r is the characteristic time scale of quenching, a scenario that supports the Kibble- 
Zurek mechanism. Interestingly, in the case of anisotropic quenching, there exists an additional 
non-adiabatic transition, in comparison to the transverse quenching case, with the corresponding 
probability peaking at an incommensurate value of the wave vector. In the special case in which 
the system passes through a multi-critical point, the defect density is found to vary as l/r 1//6 . The 
von Neumann entropy of the final state is shown to maximize at a quenching rate around which the 
ordering of the final state changes from antiferromagnetic to ferromagnetic. 

PACS numbers: 73.43.Nq, 05.70.Jk, 75.10.Jm 



I. INTRODUCTION 

A quantum phase transition corresponds to a funda- 
mental change in the symmetry of the ground state of a 
quantum system when the strength of the quantum fluc- 
tuations is appropriately tuned at zero temperature 1 -^. 
In a quantum system, statics and dynamics are inter- 
mingled, and a quantum critical point is therefore as- 
sociated with a diverging correlation length as well as a 
diverging relaxation time. The diverging time scale plays 
a non-trivial role when the system is driven through the 
quantum critical point at a uniform rate 3 -^^!^. This im- 
plies that no matter how slow the quenching may be, the 
dynamics of the system fails to be completely adiabatic 
when a quantum critical point is crossed. The possibility 
of experimental studies of non-equilibrium strongly cor- 
related quantum systems 8 - has paved the way for rigor- 
ous theoretical investigation s 9 ' 10 ! 11 ! 12 ! 13 of the dynamics 
of various model Hamiltonians when swept through their 
quantum critical points. Recently a general analysis has 
been carried out of the effects of quenching in gapless 
systems^. 

The above mentioned non-adiabaticity arising due to 
the diverging relaxation time near the critical point leads 
to the production of topological defects. The spatial dis- 
tribution of the spins in the final state is quite complex 
in comparison to the situation in which the dynamics is 
adiabatic for the entire range of time. The rate of pro- 
duction of defects can be quantified using the prediction 
of the Kibble-Zurek (KZ) theory extended to quantum 
spin chain s 15 ! 16 . If the gap of a one-dimensional quan- 
tum Hamiltonian is changed linearly as i/r, where r is 
the characteristic time scale of quenching, the dynamics 
is adiabatic for almost the entire span of time, except 



for a region in the vicinity of the quantum critical point 
called the "impulse" region where non-adiabatic transi- 
tions dominate^. The system enters the impulse region 
at a time t when the rate of change of the gap is of the 
order of the relaxation time of the system. One can show 
that t ~ yfr which eventually results in the density of 
defects decreasing as 1/y/r. In other words, there is, on 
average, a single defect in a region of length £ which also 
scales as y/r. A smaller quenching rate (larger value of 
r) yields a larger £ and the defect production is less. A 
detailed analysis of the above mechanism for exactly solv- 
able quantum spin chains is presented in Refs. 3 .4.5.6 .7 
and[l7|. Recently, the above studies have been extended 
to explore the dynamics of disordered transverse Ising 
chains 18,19 . The quenching behavior of Bose-Hubbard 
models has also been explored in recent year s 20 ' 21 . 

There is a recent upsurge in the study of quantum dy- 
namics from the point of view of an optimization prob- 
lem where the strength of the quantum fluctuations is 
quenched from a very high value to zero in order to ar- 
rive at the true ground state of a frustrated classical sys- 
tem. This approach of adiabatic quantum computation 
is popularly known as "quantum annealing" 22 ' 23 ' 24 . The 
measure of non-adiabaticity in this approach is given by 
the residual energy E TCS = Ea n — E c \, where is the 
energy of the final state, and E c \ is the energy of the true 
classical ground state. The two measures of the degree 
of non-adiabaticity, the density of defects in the previ- 
ous approach and the residual energy in the quantum 
annealing approach, are proportional to each other for a 
disorder free system. In the present literature on quan- 
tum dynamics, the terms "quenching" and "annealing" 
are often used synonymously. 

In this work, we study the dynamics of an anisotropic 
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transverse XY spin-1/2 chai n 25 ' 26 ' 2 which is driven 
across various quantum critical lines at a steady and 
finite rate. In the transverse quenching scheme, the 
transverse field is varied from — oo to oo^^ 7 -, whereas in 
the "anisotropic quenching" , which constitutes the main 
theme of this paper, the interaction term is quenched 
from — oo to oo keeping the transverse field unchanged. 
We study the "anisotropic" quenching scheme in detail, 
and compare the results with those of the transverse 
quenching scheme. In both cases, the dynamics is ex- 
actly solved via a mapping to an equivalent Landau- 
Zener problem 2 ^ through a Jordan- Wigncr transforma- 
tion from spins to fcrmions 25 . 

The paper is organized in the following way. In Sec- 
tion II, we discuss the Hamiltonian and the corresponding 
phase diagram. For the sake of completeness and to fix 
our notations and terminology at the outset, we include 
a brief discussion of the Jordan- Wigner transformation 
and the diagonalization of the spin Hamiltonian in the 
fermionic representation. We report the results of the 
transverse as well as the anisotropic quenching scheme in 
Section III. In Section IV, we discuss the behavior of the 
von Neumann entropy and the magnetization of the final 
state as a function of the quenching time. 



II. THE MODEL AND THE PHASE DIAGRAM 

The Hamiltonian of the one-dimensional anisotropic 
spin-1/2 XY chain in a transverse field is given b y 25 ' 26 ' 27 

H = ~ \ ( J *«+i + J v«+i + h <), (!) 



where the cr's are Pauli spin matrices satisfying the usual 
commutation relations. The strength of the transverse 
field is denoted by h, and J x — J y is the measure of the 
anisotropy of interactions in the x and y directions. In 
this work, J x , J y and h are chosen to be non-random. To 
explore the excitation spectrum of the Hamiltonian, we 
choose time-independent values of the parameters in this 
introductory section. In the limit J y = 0, Hamiltonian 
in Eq. (Q]) reduces to the transverse Ising Model 2 - 9 -, while 
for J x — J y it describes an isotropic XY model 2 ^. 

The Hamiltonian in Eq. |T]) can be exactly diago- 
nalized using the Jordan- Wigner transformation which 
maps a system of spin-l/2's to a system of spinless 
fermion o 25 ' 27 ' 30 . The Jordan- Wigner transformation of 
spins to fermions is given by 




= i i i ^ ] I- 1 )'* <?n, (2) 

K j=-co 

where = (°n ^ *°n)/2 are the spin raising (lower- 
ing) operators. The operator er^ is expressed in terms of 



fermion operators as cr* = 2cJ i c n — 1; thus the presence of 
a fermion at site n corresponds to a spin-up state. In the 
fermionic language, the above Hamiltonian can be rewrit- 
ten in Fourier space with a periodic boundary condition 
as 

H = - { [{J x + J y ) cos k + h] {c\c k + c t _ / .c_ fc ) 



k>0 



+ i(J x - J y ) sin k (4 c L fc - C-fcCfc}. (3) 



This Hamiltonian is quadratic in the c operators and can 
therefore be diagonalized using the standard Bogoliubov 
transformation; we then arrive at an expression for the 
gap in the excitation spectrum given b y 25 ' 27 

e k = [h 2 + Jl + Jl + 2h{J x + Jy)cosk+2J x J y cos2k] 1/2 .{A) 

The gap given in Eq. ([4} vanishes at h = T(Jx + Jy) for 
wave vectors k — and 7r respectively; it turns out that 
this signals a quantum phase transition from a quantum 
paramagnetic phase to a ferromagnetically ordered phase 
in which a discrete Z2 symmetry of the Hamiltonian in 
Eq. (al -> al -> -a y n and a* -> u z n for all n) 

is spontaneously broken. This transition belongs to the 
universality class of the transverse Ising mode l 27 ' 29 and 
is therefore referred to as the "Ising" transition. The 
spectrum is also gapless in the limit when the anisotropy 
vanishes, J x — ► J y , provided that \h/2J x \ < 1. The 
line J x = J y marks the phase boundary between two 
ferromagnetically ordered phases denoted by FM X and 
FM y , and the corresponding phase transition is called 
the "anisotropic transition" . In the FM X phase, J x > J y 
and hence the ferromagnetic ordering is in the x direc- 
tion, while it is the other way around in the FM y phase. 
One can also check that the long-range order in the FM X 
or FM y phase only exists for a relatively weak transverse 
field lying in the range —J x — J y < h < J x + J y . Each 
ferromagnetically ordered phase is further divided into a 
commensurate and an incommensurate region with the 
incommensurate wave vector fcg given by 



cos fcn = — 



h(J x + Jy) 
4(/ xJy 



(5) 



On the anisotropic phase boundary J x — J y , the in- 
commensurate wave vector therefore has a value fco = 
cos _1 (— h/2J x ). The boundary between these two re- 
gions inside a ferromagnetic phase is given by the relation 
h/(J x + Jy) = ±(l-7 2 ), with 7 ee (J x -Jy)/{J X + J y ), as 
shown by the thick dashed lines in Fig. 1. The exponents 
associated with the anisotropic transition are different 
from the Ising case and are identical to the exponents of 
a pair of decoupled Ising model o 26 ' 27 . 

We note that the points corresponding to J x — J y and 
h = ±2J y are multi-critical points since more than one 
phase boundary passes through them. We will see later 
that the quenching dynamics shows a different behavior 
when the system passes through those points. 
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FIG. 1: The phase diagram of the anisotropic XY model in a 
transverse field in the h/(J x + J y ) — 7 plane. The vertical 
bold lines given by h/(J x + J y ) — ±1 denote the Ising tran- 
sitions. The system is also gapless on the horizontal bold line 
Jx — J y for \h\ < Jx + Jy FM X (FM y ) is a long-range ordered 
phase with ferromagnetic ordering in the x (y) direction. The 
thick dashed line marks the boundary between the commen- 
surate and incommensurate ferromagnetic phase. The thin 
dotted lines indicate the adiabatic and impulse regions when 
the field h is quenched from —00 to 00. 



III. QUENCHING SCHEME AND RESULTS 

In this section, we will discuss the results obtained 
for the quenching dynamics of the Hamiltonian in Eq. 
(fl]) for two different schemes. In the first scheme, the 
time-dependent transverse field is of the form h(t) = 
t/r, where t is varied from —00 to 00, and r is the 
characteristic time scale of quenching, often referred to 
as the "quenching time" . This "transverse quenching" 
scheme has been studied extensively for the XY chain by 
Cherng and Levito^, and for the transverse Ising case 
by Dziarmaga4 and by Polkovnikov 7 . We shall briefly 
present the results for the transverse quenching which 
will be helpful in comparing it with the results obtained 
in the other scheme, namely, the "anisotropic quench- 
ing" . 

As mentioned in the Introduction, in the "anisotropic 
quenching" scheme, the interaction term J x (— t/r) is 
quenched from a very large negative initial value to a 
very large positive final value, with J y and h held fixed at 
some positive values. For our numerical studies, we will 
set 2 J y > h for the reason explained below. In the initial 
state, all the spins are antiferromagnetically ordered in 
the x direction. On the other hand, the final state should 
correspond to a state with a perfect ferromagnetic order 
in the x direction had the dynamical evolution been adi- 



abatic for the entire span of time. However, due to the 
non-adiabatic transitions near the quantum critical re- 
gions, the final state in the limit t — ► 00 will not be per- 
fectly ordered and will include a finite fraction of kinks 
or spins aligned anti-parallel to the direction of ordering. 

The Hamiltonian in Eq. ([3]) decouples into a sum of in- 
dependent terms, H(t) = ^2 k>0 Hk(t), where each Hk{t) 
operates on a four-dimensional Hilbert space spanned 
by the basis vectors |0), \k) = c[\0), \ - k) = cL fc |0), 



and \k,—k) 



t„t 



cic 



|0). The vacuum state where no c- 



particle is present is denoted by |0) which corresponds 
to a spin configuration with all spins pointing in the 
—z direction. The occurrence of only bilinear terms like 
c l c -fe m t ne Hamiltonian in ([3]) ensures that the parity 
(even or odd) of the total number of fermions given by 
n k = c\ck + c_ k C-k is conserved for each value of k > 0. 
Thus the states |0) and \k, — k) are coupled to each other 
by the Hamiltonian, while the states and | — k) remain 
invariant. 

To study the dynamics of transverse quenching, it is 
sufficient to project the Hamiltonian H k (t) to the two- 
dimensional subspace spanned by |0) and \k,—k) since 
the ground state of Eq. (J5]) for each value of k lies within 
this subspace. In this subspace, the Hamiltonian takes 
the form 



H k (t) 



[h + (J x + J y ) cos k] I 2 

(J x + Jy) COS k 



i(Jx 

i(J x —Jy)smk —h — (J x + Jy) cos k 



- J y ) sin k 



where I2 denotes the 2x2 identity matrix. A state in 
this subspace can be represented as a linear superposition 
ipk{t) = Ufc(i)|0 > +Vk(t)\k, — k >, where the amplitudes 
Uk(t) and Vk{t) are time-dependent. The initial condition 
in the transverse quenching scheme is given by itfc(— 00) = 
1 and v k (— 00) = 0. The time evolution of a generic state 
is governed by the Schrodinger equation 



idtMt) = H k {t) Mt)- 



(6) 



The projection of the Hamiltonian to the 2x2 Hilbert 
space has effectively reduced the many-body problem to 
the problem of a two-level system. The off-diagonal term 
of the projected Hamiltonian, A = (J x — J a )sinfc, rep- 
resents the interaction between the two time-dependent 
levels Ex,2 — + (Jx + Jy) cosfc]. The Schrodinger 

equation given above is identical to the Landau- Zencr 
problem of a two-level system, where the off-diagonal 
terms of the Hamiltonian determines the non-adiabatic 
transition probability pk for the wave vector k. Using the 
Landau-Zener transition formula, this transition proba- 
bility of excitations at the final time is given by2^ 



Pk 



-27T7 



(J) 



where 7 = A 2 / |t§(-Ei — E 2 )\. Equivalently, p k deter- 
mines the probability that the system remains in the ini- 
tial state |0) at the final time. For the transverse quench- 
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ing case, one can use the relation in Eq. Q to obtain p k 
as a function of the anisotropy term J x — J y and r as£ 



p k = e-^iJx-Jy) 



(8) 



This leads to an expression for the density of kinks n 
generated due to non-adiabatic transitions, 



Pk 



K*jT \J X — Jy 



(9) 



Eq. §§§ shows that the kink density decreases as l/\ff 
for large r as predicted by the Kibble-Zurek theory and 
proved by Cherng and Levitov in the case of transverse 
quenching^. 

We shall now focus on the anisotropic quenching 
scheme where the interaction term J x {t) = t/r is changed 
from -co to oo, with J y and h held fixed at some posi- 
tive values. Let us first point out the range in time where 
the system fails to follow the instantaneous ground state 
(namely, the region where the relaxation time is large or 
divergent) as J x is varied. If J x {t\) + J y = —h, i.e., for 
t\ = —r(h + J y ), the system undergoes an Ising transi- 
tion from the initial antiferromagnetic phase to a para- 
magnetic phase. When J x is further increased so that 
Jxit^) + J y — h, i.e., for t 2 = (h — J v )t, there is a phase 
transition from the paramagnetic to FM y phase. It is to 
be noted however, that for J x — — J y , there is no fur- 
ther anisotropic transition since the magnitude of h is 
greater than J x + J y = 0, and the system stays param- 
agnetic. Eventually, at a time £3 given by J x {t^) = J y , 
with h < 2J y , the system evolves to the FM X phase. The 
non-adiabaticity dominates in the vicinity of these three 
quantum critical points, i.e., for t close to tx,t 2 and £3. 
One should also note that if we fix h > 2 J y at the outset, 
there cannot be any anisotropic transition, and the sys- 
tem directly evolves from the paramagnetic phase to the 
FM X phase through the Ising transition only; hence there 
are two, rather than three, regions of non-adiabaticity. 
All the transition points are depicted in the static phase 
diagram of the model shown in Fig. 1. The adiabatic 
and non-adiabatic regions of time evolution are shown 
schematically in Fig. 2. 

As mentioned already, in the anisotropic case, the re- 
duced Hamiltonian in the |0), \k, — k) subspace includes 
a time-dependent J x term with static positive values of 
J y and h. The eigenstates of the Hamiltonian in the limit 
t — * ±00 are 



and 



|e lfc ) = sin(fc/2)|0) + £cos(fc/2)|fc, -k) 



|e 2fc ) = cos(fc/2)|0) - isin(k/2)|Jfc, -k), 



with eigenvalues Ai = t/r and X 2 = — t/r respectively; 
the system is in the state \e\ k ) initially. A general state 
vector can be expressed as a linear combination of \e)\ k 
and |e) 2 fc, 



FMy-> FM, 
Impulse 




-(h + J y )T 



(h-J y )x 



Time 



FIG. 2: Schematic diagram showing the divergence of the re- 
laxation time of the quantum Hamiltonian at the quantum 
critical points. As discussed in the text, the dotted verti- 
cal lines denote the impulse region of dynamics where non- 
adiabatic transitions play a prominent role. 



The initial condition in the anisotropic case is 
C\k{— 00) = 1 and C2fe(— 00) = 0. The amplitudes Uk 
and Vk are related to the new coefficients as 



Uk = Cife sin(fc/2) + C 2 k cos(fc/2), 
v k = iCik cos(fc/2) - iC 2 k sin(fc/2). 



(11) 



The Landau-Zener transition probability can be obtained 
by numerically solving the Schrodinger equation in Eq. 
^ using the basis vectors |0) and \k,—k), with the appro- 
priate initial conditions for u k and v k - The non-adiabatic 
transition probability p k at the final time is simply given 
by Pk = |Cife(oo)| 2 . We have shown the variation of p k 
with k obtained through numerical integration in Fig. 3. 

The situation is apparently complicated in the 
anisotropic quenching case because the off-diagonal terms 
of the 2x2 Hamiltonian matrix in the |0) and |fc, — k) 
basis are time-dependent. However, by an unitary trans- 
formation to the new set of basis vectors \e)\k and \e) 2 k, 
one can get rid of the time-dependence of the off-diagonal 
terms and map the anisotropic quenching problem to an 
equivalent Landau-Zener problem. The unitary transfor- 
mation is given by H' k (t) — W Hk(t)U , where 



U = 



cosO/2) sin(fc/2) 
-isin(fc/2) icos(fe/2) 



\Mt)) =C lk (t)\e lk )+C 2k {t)\e 2k ). 



(10) 



and the new Hamiltonian is 
H' k (t) = - [h + {J X + Jy) cos k]I 2 

Jx + Jy cos 2k + h cos k J y sin 2k + h sin k 
J y sin 2k + h sin k —J x — J y cos 2k — h cos k 

Therefore, the anisotropic quenching is now placed on 
the same footing as the transverse quenching case, with 
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FIG. 3: pk vs k as obtained numerically and analytically for r 
equal to 0.5 and 5. We have fixed J y = 1 and h — 0.2. In the 
region marked r = 0.5, the solid line is numerical and dashed 
line is analytical. In the region marked r = 5, the dashed line 
with smaller spacing is numerical and dashed line with larger 
spacing is analytical. 



off-diagonal terms in the Hamiltonian matrix which are 
time-independent. The time evolution of the amplitudes 
C\k{t) and C2fc(i) is dictated by the Schrodinger equa- 
tion, 



.dC u 
i — — 
dt 



dG 



2k 



dt 



(J x (t) + J y cos 2k + h cos k) C\h (t) 

+ (J y sin 2k + h sin k) Cik (t) 

(J y sin 2k + h sin k) C\k (t) 

- (Jx(t) + Jy cos2fc + hcosk) C 2 fc(t), (12) 



where we have dropped the effect of the term involving 
the identity matrix in the Hamiltonian since this affects 
Cifc and C2k by the same time-dependent phase factor, 
and will therefore have no effect on the density matrix. 
The non-adiabatic transition probability now depends on 
J y and h, and is given by 



Pk = |Cifc(oo)| s 



g — 7rr(J y sin 2k-\-h sin k) 



(13) 



In Fig. 3, we plot -p k as a function of the wave vector 
k along with the numerical results. We see that in con- 
trast to the transverse case, there is an additional region 
of non-adiabaticity peaked at the incommensurate wave 
vector fco = cos~ 1 (—h/2J y ). This non-adiabaticity arises 
due to the existence of the anisotropic transition of the 
underlying static XY model at J x (t) = J y . The situ- 
ation with zero transverse field is a special case where 
the regions of non-adiabaticity peak at k — 0, tt/2 and 7r, 
respectively. 

As in the transverse quenching case, we now measure 



the density of kinks given by 



w dk 



Pk 



71 dk 



-7tt( J y sin 2k-\-h sin ky 



■ (14) 



In Fig. 4, the plot of the kink density as a function of 
the quenching time r is shown, which clearly shows that 
n oc \j\fr for large values of r. This finding supports 
the prediction of the Kibble-Zurek mechanism even in 
the anisotropic quenching case. One can also find an 
approximate analytical form of n in the following way: for 
large t, only the modes very close to the critical modes 
contribute. For k — » and tt, pk can be approximated 
as exp[— ttt(2J v + h) 2 k 2 ] and exp[— ttt(2J v — h) 2 (-K — k) 2 ] 
respectively. With this p k , the density of kinks produced 
by the modes near k — and tt is given by 



2tt\/t 



1 



1 



2J y - h 



2 1/71 



7T^?(4J 2 -h 2 )- 



By expanding around fco, we find that the contribution 
to n from modes with k ~ ko is equal to the contribution 
from k = and k — tt taken together, 



n-2 



* dk 



1 



-7tt (2 J y cos 2kQ-\-h cos ko) 2 (k—ko) 2 



1 



TTy/r \2Jy cos 2&o ~t~ h cos fco I 
Therefore the total kink density is given by 



m + n-2 



4J„ 



(U 2 -h 2 \ 



(15) 



Fig. 4 shows that this approximate form embraces the 
exact result perfectly in the limit of large r. 




FIG. 4: Variation of kink density n with r as obtained by 
numerical integration of Eq. (fTI]) for h = 0.2,0.4,0.6,0.8 
(from bottom to top in the large r region), with J v = 1. For 
large r, n increases with increasing h, whereas for small r, it 
decreases with increasing h as shown in the inset. 
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Using Eq. (fT5|) . one can find the variation of the kink 
density with h, for large values of r, as 



dn 
dh 



8hJ.„ 



'K\Jt (4J2 - h 2 ) 2 • 



The slope is positive for all values h indicating an in- 
crement in kink density with increasing h as shown in 
Fig. 4. On the other hand, for small r, the density 
n is found to decreases with increasing h. In the limit 
7rr( Jy sin 2/c + ft, sin A:) 2 <C 1, one can expand the ex- 
ponential in pk in Eq. (1141) . retaining only terms up 
to first order in r. One can then show that dn/dh 
is negative. The crossover from the small r to the 
large r behavior occurs around a typical quenching time 
r = ln2/[7r(J„ sin2fc + hsmk) 2 }. The significance of the 
crossover time is explained below. 

To derive a characteristic time scale of the anisotropic 
quenching, we first note that the minima of the proba- 
bility pk occurs at a wave vector value k, where 



-h±Jh 



cos A; = 



32J 2 



8 J v 



(16) 



The presence of a non-zero transverse field h leads to an 
asymmetry in pk on the either side of the maxima at ko 
as is evident at smaller values of r (Fig. 3). We now 
define a time scale To so that at r = To, pk = 1/2 at the 
minima 6 k — k. This implies 



In 2 



to 



ir{J y sin 2k + h sin k) 2 



(17) 



The two conjugate values of k given in Eq. ([16]) lead to 
a pair of ro's for the anisotropic quenching, in contrast 
to the transverse quenching case. For h — 0, these two 
values coalesce into one. The crossover of density shown 
in Fig. 4, takes place roughly around the smaller value 
of To denoted by To2- For t 3> toi,to2, the dynamics is 
nearly adiabatic except for the modes very close to the 
critical modes (Fig. 5). 

Eq. (Ti"5|) is valid for 2J y > h. On the other hand, if 
we take h > 2J y , the system evolves directly from the 
paramagnetic to the FM X phase; there is no anisotropic 
transition. In this case, Eq. (fl"4"|) gives the kink density 
to be 



1 



2tt^ 



1 



1 



h + 2J v h- 2J V 



7T^ (ft 2 "4 J 2 )' 



(Note that dn/dh is negative in this case). The trans- 
verse Ising model result given in Ref. |4j can be obtained 
by setting J y = in the above expression. 

Finally, we consider what happens if 1 — 2J y /h is close 
to zero. If 2J y = ft, the system passes through a multi- 
critical point when J x (t) equals J y ; we therefore expect 
an unusual behavior of the kink density n in this case. If 
2 J y = ft, and r is large, we find that Eq. (fl4|) is dominated 
by the region near k = 7r; since pk can be approximated 
there by exp[— 7tt/i 2 (7t — fc) 6 /4], the contribution of this 




FIG. 5: Variation of pk vs k for r = 5, toi,to2 and 0.01 from 
bottom to top, where pt = 0.5 at the two minima for r = roi 
and To2- We have fixed J v = 1 and h = 0.2. For r S> roi,To2, 
the evolution is almost adiabatic except at k — 0, ir and ko, 
where the system remains frozen in its initial state. For very 
small t, pk is close to 1 for all k. 



region to n goes as 1/(/i 1 / 3 t 1 / 6 ). [The contribution from 
the region near k = to n goes as l/(/iy / r) which is 



much smaller than 1/(/i 1 / 3 t 1 / 6 ) 



if h 2 r > 1.] Now, if 



1 — 2J y /h is non-zero but small, pk can be approximated 
near k = 7T by exp[-7rr{(/i - 2 J y ) (tt - k) + h(n- fc) 3 /2} 2 ] . 
We then see that the term T(h—2J y ) 2 ('K—k) 2 will become 
more important than the term rh 2 {it — k) 6 and the kink 
density will show a crossover from a 1/t 1 / 6 behavior to 
a 1/s/r behavior when r increases beyond a value which 
is of order h/\h— 2J y | 3 . 



IV. ENTROPY AND MAGNETIZATION OF 
THE FINAL STATE 

The quenching dynamics of the spin model is dictated 
by the Schrodinger equation and is obviously unitary. 
Therefore, the final state must be a pure state described 
by a density matrix of the product form p = Pk, where 
p k is given by 



Pk 
Ik 



qk 
1 - Pk 



(18) 



and qk — Cifc(oo)C2 fe (oo). The diagonal elements of the 
reduced 2x2 density matrix are smooth functions of 
k and are independent of the total time of quenching 
(which we now set equal to T for convenience, with T 3> 
1), whereas the off-diagonal terms are rapidly oscillating 
functions of both k and T. 

The final state, even though a pure state, has a fairly 
complicated local structure; its local properties in the 
real space is identical to that of a mixed state with a 
finite entropy 6 -. In other words, the final state can be 
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viewed as a superposition of different configurations of 
magnetically ordered domains. The off-diagonal terms of 
the density matrix pk can be made to vanish at the final 
time T upon coarse-graining in the wave vector k. The 
final state therefore may be viewed locally as a mixed 
state described by a decohered reduced density matrix 
p D given by 



To quantify the amount of information lost in the deco- 
herence process, we consider the von Neumann entropy 
density of the system, s = — tr(p£> In pd), namely, 

f" dk 

s = ~ — [PkHPk) + (l-pfe)ln(l-p fc ) ]. (20) 
Jo v 

In the case of the transverse quenching 6 , the maximum 
of the entropy density occurs near the value of r = 2r 
where the minimum value of pu = 1/2. In the anisotropic 
case, as mentioned already, tq is non- unique (Fig. 5). 
The entropy density s when plotted against t/to2 shows 
a maxima at a quenching rate r ~ 2to2 (see Fig. 6), while 
no special behavior is noted near tqi] this establishes T02 
as the characteristic time scale of the anisotropic quench- 
ing. For sufficiently small r, non-adiabaticity dominates 
and the system stays in its initial state with high prob- 
ability, namely, the system more or less retains its ini- 
tial antiferromagnetic ordering even at the final time. 
The system looks like a pure state even at small length 
scales. In the large r limit, on the other hand, the sys- 
tem evolves essentially in an adiabatic way to the final 
ground ground state, and the non-adiabatic transition 
probability is small. The final state is ferromagnetically 
ordered. Therefore, the entropy density vanishes asymp- 
totically for both large and small values of r, and the 
maximum occurs at an intermediate value of r ~ 2to2 
where the defect in the local ordering of the final state is 
maximum. 

The above arguments can be justified further by ex- 
ploring the magnetization of the system in the final state. 
Let us recall that the initial state is antiferromagnetically 
oriented in the x direction. The average energy of the 
mode k at the final time is related to pk as (2pk — 1) J x ] 
hence, integrating over all the modes gives the total en- 
ergy per site to be (2n — 1) J x , where n is the density 
of kinks defined above. On the other hand, at the final 
time, J x 3> J y and h; hence the effective Hamiltonian 
is H = — J2 n i^K+i' an d the total energy density is 
— J x rr? x . Putting all this together and using Eq. (fT5|) , we 
arrive at an expression for the magnetization, 

/ g J \ V 2 

"^"K 1 " ^^-*) ) ■ (21) 

and m x = whenever the right hand side of Eq. (f2Tj) is 
imaginary. In Fig. 6, we plot the final magnetization as 
a function of t/tq2 which shows that the magnetization 
starts to become non-zero when the quenching time is of 



the order of 2to2- For slower quenching (larger values of 
t), the dynamics tends to be more and more adiabatic, 
and the final magnetization monotonically increases to- 
wards the saturation value of unity. 

Let us now shift our attention to the behavior of the 
staggered magnetization m sx . In the limit of small r, 
the system fails to follow the instantaneous Hamilto- 
nian; hence the final state retains an antiferromagnetic 
order with vanishing total magnetization. In other words, 
the total magnetization of the even sites cancels the to- 
tal magnetization of the odd sites. Using similar argu- 
ments as given above, the staggered magnetization can 
be shown to behave as (2n — l) 1 / 2 ; hence it vanishes at 
n = 1/2 and stays at zero for smaller values of n (i.e., 
higher values of r). On the other hand, for smaller val- 
ues of r, n > 1/2 and there is a non-zero value of m sx . 
Fig. 6 shows that the staggered magnetization of the fi- 
nal state vanishes at a quenching rate close to 2tq where 
the ferromagnetic order starts to set in. 

We can use similar arguments as given in the case of 
entropy density to understand the variation of the mag- 
netization with t as presented above. We have argued 
already that even though the final state is a pure state, 
it can also be viewed locally, or with a coarse-grained 
wave vector scale, as a decohered (mixed) state. If the 
quenching time scale t <C T02, there are very few local 
defects and the system retains the initial antiferromag- 
netic order. On the other hand, for r 3> tq2, the final 
state is locally ferromagnetically ordered. Hence, there 
exists an intermediate region of r ~ T02 where the an- 
tiferromagnetic order decreases rapidly. Eventually, for 
r > 2to2, the final state starts to acquire ferromagnetic 
ordering. 



V. CONCLUSIONS 

We have studied the adiabatic quantum dynamics of 
an anisotropic XY spin-1/2 chain in a transverse field 
when the parameters of the Hamiltonian are quenched 
at a steady and finite rate. Our emphasis lies on a new 
scheme, namely, the anisotropic quenching, where the 
interaction strength in the x direction is quenched from 
—00 to 00 at a uniform rate dictated by a time scale r. 
At first sight, the Landau-Zener theory does not seem to 
be applicable in this case due to the presence of time- 
dependent off-diagonal terms in the reduced 2x2 Hamil- 
tonian. However, through a unitary transformation to 
an appropriate basis, the time dependence can be en- 
tirely shifted to the diagonal terms, and Landau-Zcncr 
theory turns out to be applicable once again. 

In the process of the anisotropic quenching, the system 
is swept across all the three quantum critical lines in 
the phase diagram if h < 2J y , and across two quantum 
critical lines if h > 2J y . In both these cases, the density 
of defects is found to vary as 1/\/t. In the special case of 
h = 2J y , the system is swept across a multi-critical point, 
and the density of defects is found to vary as 1/r 1 / 6 . 
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FIG. 6: Variation of von Neumann entropy density s, stag- 
gered magnetization m sx and magnetization m x as a function 
of t/tq2, for J y — 1 and h — 0.2. 



Finally, we find a characteristic quenching time T02 
around which the von Neumann entropy of the final state 
maximizes and ferromagnetic ordering starts to set in. 
The residual energy of the final state is proportional to 
the defect density and hence scales as 1 / y/r. 
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